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^  Mathematical  formulae  are  given  Cron  which  target  strength  spectra 
and  acoustic  intensity  vectors  are  calculated.  The  excitation  is  a  time— 
harmonic  plane  wave,  intensity  vector  plots  show  the  considerable  distortion 
of  the  incident  sound  field  that  may  be  present  when  an  elastic  scatterer  is 
excited  at  a  resonant  frequency. 
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Itw  dynamics  and  acoustics  of  a  uniform  spherical  Shall  is  a  basic 
problem  in  the  field  of  fluid-structure  Interaction  because  the  relative 
simplicity  of  the  geometry  and  equations  enables  closed-form  solutions  in 
tecw  of  series  of  Legendre  polynomials  and  spherical  Bessel  functions,  thus 
rendering  possible  general  Insight  into  the  effect  of  fluid  loading  on -the 
dynamics  of  non-planar  elastic  surfaces.  Additionally,  numerical  results  maq 
be  used  as  a  check  of  the  correctness  of  computer  programs  based  on  purely 
numerical  formulations  such  as  the  coupled  finite  element  and  Halaholts 
integral  equation  methods.  Much  of  the  theoretical  work  on  the  dynamics  and 
acoustics  of  a  thin  spherical  shell  is  contained  in  the  standard  text  [1]  on 
fluid-structure  interaction.  Nong  a  Hayek  [2]  have  analysed,  both  analyti¬ 
cally  and  experimentally,  the  vibration  and  acoustic  radiation  from  point 
force  excited  shells  in  air  and  in  water.  They  found  good  agreement  between 
the  resonant  frequencies  determined  by  theory  and  experiment.  Other 
references  of  interest  are  contained  in  the  locations  cited. 


—  >>  Acoustic  intensity  vector  plots  can  illustrate  vividly  the  power  flow  in 
an  acoustic  fluid.  \  Spicer's  plots  [3]  show  the  subsonic  and  'leaky*  waves 
associated  with  the  interaction  between  an  elastic  plate  and  a  fluid  when  the 
former  is  driven  toy  a  time-harmonic  force.  James  [4]  gives  vector  plots  Which 
Show  that  there/Hty  be  significant  interchange  of  power  between  a  duct’s  wall 
and  its  contained  fluid  Mien  the  excitation  is  a  force.  Intensity  vector 
plots  of  thSscattered  sound  field  due  to  plane  wave  excitation  of  an  infinite 
fluid-filled  cylindrical  shall  have  been  given  by  janes  [5].  vector  plots  of 
the  total  sound  field  were  also  obtained  but  not  reported.  The  plots  did  not 
Show  Significant  distortion  of  the  incident  wave,  except  for  the  special  case 
of  S  pressure  release  scatterer,  possibly  because  there  was  no  marked  resonant 
response  of  the  shell. 


— '  It  Is  the  main  purpose  of  this  memorandum  to  examine,  by  means  of 
acoustic  intensity  vector  plots,  the  distortion  of  the  Incident  sound  field 
due  to  the  presence  of  a  spherical  Shell  Muse  target  strength  spectrum 
contains  well-defined  resonant  peaks. The  mathematics  is  given  in  Section  2 
with  a  minimum  of  explanation  because  most  of  it  can  be  found  scattered 


throughout  the  standard  text  [1]. 
obtained  from  Fortran  programs  are 


-strength  spectra  and  vector  plots 
sad  in  Section  3. 
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The  radial  and  tangential  displacements  at  ths  shell's  surface  are 
seed  as  a  series  of  Legsndre  polynomials  of  argtmmnt  n  -  cose,  vis. 
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In  tha  abova,  h  is  tha  shall' a  thicknass  and  a  1s  its  naan  radius; 

E i -Eh/ [ ( l-v2 )*2 ]  whara  E  is  Young's  aodulus  and  v  is  Poisson's  ratio; 

0*  -  h2/12a2  and  Xn  *»  n(  n+1 ) ;  ps  is  tha  dansity  of  tha  shall' s  aatarlal.and 
P 1  and  p2  ara  tha  dansitlas  of  tha  intarior  and  axtarior  fluids  whoaa  sound 
valodtias  ara  ct  and  c2,  raspactivaly;  kl  and  k2  ara  acoustic  wacva  nunbara, 
u/Cj  and  m/c2;  j„  and  hn  ara  spharlcal  Bassal  functions,  and  Pn  ara  Isgandra 
polynaaials  of  dagraa  n;  tha  prias  on  various  quantltias  danotas  diffbrantla 
tion  with  raspact  to  thair  argiaaants.  Tha  tina  harmonic  factor,  axp( -ii.it), 
will  ba  caalttad  throughout.  Tha  gaaastry  of  tha  problan  is  shown  in  Figuxa  l 
in  which  it  Is  avldant  that  tha  problan  is  axis ynwa trie,  i.a.  in  a  syatan  of 
spharlcal  ooordlnatas  (R,e,$)  tha  displacaannts  ara  lndapandant  of  tha 
drcunfbrantial  angla  a. 


In  aquation  (3)  Fn  ara  tha  coafficiants  of  tha  Isgandra  sarias  axpanalon 
of  tha  amcitatlon  strass 
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n-0 


Is  the  acoustic  pressure  due  to  the  Incident  plane  wave  and 
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Is  the  pressure  scattered  by  the  shell  behaving  as  a  rigid  laanvable  body. 

Coafcining  equations  (7)  and  (8),  with  the  help  of  the  Wronsklan  relation 
jn(x)y'n(x)<-j,n(x)yn<x)  -  1/x2,  gives  the  excitation  as 
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the  pressure  in  the  exterior  fluid  Is  the  sub  of  three  tens,  vis 

* 

p2<R,0)  -  pt(R,0)  +  Pr(R,0)  +  P«( R, 8 )  (10) 


the  first  ten,  p*,  Is  the  pressure  due  to  the  incident  wave,  equation 
(7)>  the  second  ten,  pr,  is  the  pressure  scattered  by  the  shell  behaving  as 
rigid  body,  equation  (8);  and  the  third  ten,  p#,  is  the  sound  field  arising 
from  elastic  vibrations  of  the  shell,  this  tan  is 
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It  is  very  easily  verified  that  the  boundary  condition 
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which  eafeodles  continuity  of  nonal  dlsplaceeent  between  exterior  fluid  and 
shell  surface  Is  satisfied,  the  effect  of  the  interior  fluid  is  coapletely 
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aaldantly  tha  radiation  loading  of  tha  axtarlor  fluid. 

Tha  scattarad  prassura  la  amply 
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which  la  put  Into  its  far-flald  fora,  pBf ,  by  zaplacing  tha  spharlcal  Baidcal 
function  by  its  valua  for  a  largo  axguaant,  vis., 
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n»  magnitude  and  direction  of  the  lntanalty  victor  In  the  XZ  plana  ara 
given  lay  the  formulae 
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(a)  gtoeral 


h  Fortran  program  has  been  written  to  calculate  tha  a  target  strength 
defined  as 


TS  -  20.0  X  log10[ ( H/Pj, ) I Ps£( R* © ) I) 


where 


Psf(R.»)  -  Ida  (Pr(R,©)  +  p#(R,©)3 


la  tha  preasure  scattered  to  tha  far- field.  Tha  intensity  vectors  ara 
calculated  from  computations  of  pressure  and  displacements  close  to  the  Shell 

In  Figures  2-8,  the  following  constants  in  si  units  were  used  to  obtain 
the  plotted  data  i 

Steel  Shellt  E  -  19.5X1010  v  -  o.29  p,  -  7700.0  a  -  1.0 
h  -  0.01/0.1  for  the  thin/thicJc  shell 


Mater i 


p  -  1000.0  C  -  1500.0 


Shall  damping  mi  included  in  the  confutations  by  netting  Young's  Modulus  to 
ths  complex  value  E< 1-iri )  where  n»  the  hyeteretic  loss-factor,  was  chosen  as 
0.01. 


(b)  Target  strength  spectra 

Figure  2  Shows  aonostatic  target  strength  spectra  of  the  thin  Shell,  the 
frequency  range  extends  to  a  ka-value  of  8.4.  In  the  low  frequency  regime 
there  are  doninant  resonant  peaks  which  are  labelled  by  their  individual 
haznonic  nuebers.  These  resonances  belong  to  the  lower  branch,  which  has  no 
resonances  at  n-0  or  1,  of  the  frequency/node-nunber  dispersion  plot  [2]. 

The  resonant  peaks  of  ths  enpty  Shell  are  particularly  high,  the  n-2  peak 
being  almost  20dB  above  the  'rigid  body*  datum  plot.  Above  400Hz  the  saooth 
variation  in  the  target  strength  spectra  is  caused  by  the  dominance  of  nodes 
which  are  excited  at  frequencies  remote  from  their  natural  frequencies.  In 
this  frequency  regime  the  plots  are  similar  to  those  of  sound  scattering  from 
an  infinite  cylindrical  shell  [S]. 

In  Figure  3,  Which  shows  target  strength  spectra  of  the  thick  shell,  the 
resonant  peaks  are  broader  and  extend  throughout  the  entire  frequency  range, 
their  levels  being  up  to  lSdB  above  the  'rigid'  datum  plot.  When  the  shell  is 
filled  with  water  there  is  a  small  peak  due  to  an  n-1  resonance,  presumably  of 
the  upper  branch  [2]  of  the  dispersion  plot. 


(c)  Intensity  Vectors 

In  Figures  4-8,  which  show  intensity  vector  plots,  in  the  Z-x  plane,  at 
selected  frequencies,  the  plotted  lengths  of  the  vectors  are  proportional  to 
I1/2.  They  have  been  normalised,  in  each  plot  separately,  to  a  maximal  length 
equal  to  the  grid  spacing.  Due  to  the  symmetry  of  the  excitation  and 
geometry,  it  is  only  the  vectors  in  the  region  X  >  0  that  have  been  plotted, 
the  plotting  region  extending  over  the  ranges  -2  <  Z  (  2  and  0  4X42. 

In  Figure  4  plots  of  a  rigid  (i— ovable)  sphere  are  shown  at  250  and 
1000HX.  There  is  a  minimum  of  distortion  of  the  incident  sound  field, 
although  the  initial  development  of  the  shadow  sons  behind  the  sphere  is 
evident  at  iooohz. 

In  Figure  5  the  large  distortion  of  the  incident  sound  field  at  resonant 
frequencies  of  the  empty  thin  shell  is  shewn.  At  2710s  ( n-2 )  there  is  a 
mas xbmuw  of  power  flow  into  the  shell  at  the  poles,  and  a  maximum  of  power  flew 
out  of  the  site  11  at  the  equator.  At  342Hz  (n-3 )  there  is  a  minimum  of  power 
flow  into  the  Shell  at  the  poles,  but  still  a  maximiwi  of  power  flow  out  of  the 
Shell  at  the  equator.  This  behaviour  is  consistent  with  the  behaviour  of  the 
Legendre  polynomials  f2  snd  P3. 

The  vector  plots  of  the  empty  thin  shell  at  788Ks  (dip  in  spectrum)  and  at 
13O0BS  (a  maximum  in  spectrum)  are  shown  in  Figure  6.  Unlike  the  previous 
plots  Which  Show  violent  distortions  at  resonance,  these  plots  are  much  the 
sane  as  ths  'rigid  body'  plot  of  Figure  4B. 

In  Figure  7,  the  vector  plot  of  the  fluid-filled  thin  Shell  at  the 
resonant  frequency  of  l94Hs  (n-2)  is  much  the  same  as  the  plot  of  Figure  SA. 


llw  vector  plot  of  tho  fluid-filled  thin  shell  at  175Hs  (dip  in  spectrum) 
Shorn  significant  distortion  at  the  equator  only. 

In  Figure  8A,  the  vectors  of  the  empty  thick  shell,  are  shown  at  497Hz, 
the  resonant  frequency  of  the  n-2  mode.  The  most  significant  feature  is  the 
the  intense  power  circulation  in  the  illuminated  region.  In  Figure  8B,  the 
vectors  of  the  fluid-filled  thick  shell  are  shown  at  556Hz,  the  resonant 
frequency  of  the  n-l  node  which  is  thought  to  belong  to  the  upper  branch  of 
the  shell  dispersion  relation.  It  is  difficult  to  Interpret  these  plots  due 
to  interference  from  neighbouring  modes. 


The  limited  number  of  intensity  vector  plots  shown  herein  have  illustrated 
vividly  the  substantial  distortion  of  the  incident  sound  field  that  may  take 
place  when  a  system  is  excited  at  a  resonant  frequency,  which  is  defined  as  a 
frequency  at  which  there  is  a  sharp  maximal  in  the  target  strength  spectra. 

Due  to  the  importance  of  the  sphere  as  a  standard  target  for  underwater 
reflectivity  measurements ,  the  following  follow-up  projects  may  be  of  sene 
values 


(i)  Development  of  Fortran  program  for  predicting  target  stength  and 
intensity  vectors  of  a  sphere  composed  of  an  arbitrary  nuaber  of 
layers  of  elastic  solids  and  fluids.  Skelton's  work  on  layered 
cylinders  [6]  is  of  relevance  here. 

(li)  Investigation  of  sound  scattering  by  a  fluid-filled  spherical 
shell  with  point  attachments  such  as  masses  and  springs.  The 
method  of  dynamic  stiffness  coupling  as  used  by,  for  example, 
Spicer  [31  enables  the  point  attachments  to  be  treated  as  forces 
of  constraint. 
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FIG.,  1  SOUNO  SCATTERING  FROM  SPHERICAL  SHELL 
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FIG. 2  MONOSTATIC  TARGET  STRENGTH  OF  THIN  SHELL 


A.  FREQUENCY  -  ZJiSz  (n-2  resonance) 
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B.  FREQUENCY  =  342  (n-3  resonance) 
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FIG. 5  INTENSITY  VECTORS  OF  THIN  EMPTY  SHELL 
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A.  FREQUENCY  -  788Hz  (dip  in  spectrum) 
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B.  FREQUENCY  *  1300Hz  (maximum  in  spectrum) 


FIG. 6  INTENSITY  VECTORS  OF  THIN  EMPTY  SHELL 
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A.  FREQUENCY  -  194Hz  (n-2  resonance) 
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B.  FREQUENCY  »  175Hz  (dip  In  spectrum) 
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FIG. 7  INTENSITY  VECTORS  OF  THIN  FLUID-FILLED  SHELL 
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B.  FREQUENCY  *  556Hz  (n=l?  resonance) 
FLUID-FILLED  SHELL 
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FIG, 8  INTENSITY  VECTORS  OF  THICK  SHELL 


